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OBJECTIVE  AND  APPROACH 


Detection  of  underground  nuclear  explosions  includes  the  spectral  analysis  of 
seismograms,  an  important  portion  of  which  is  the  contribution  of  Love  waves.  Field 
evidence  suggests  that  it  may  be  possible  to  discriminate  between  nuclear  events  and 
earthquakes  by  examining  the  Love  wave  records.  The  spectra  for  these  events  are 
different  because  an  earthquake  generates  shear  waves  directly,  whereas  an  underground 
explosion  generates  P-waves,  from  which  Love  waves  are  produced  by  scattering  from 
material  heterogeneities. 

The  objective  of  this  research  is  to  investigate  experimentally,  and  support  by 
theoretical  calculations,  the  formation  and  propagation  of  Love  waves  from  a  P-wave 
source  due  to  scattering  at  material  heterogeneities.  The  approach  is  shown  schematically 
in  Figure  1.  In  these  experiments,  a  spherical  piezoelectric  crystal  (P-wave  generator)  is 
cast  in  a  surface  layer  of  rock  simulant  overlying  a  higher  impedance  granite  substrate. 
Exitation  of  the  piezoelectric  crystal  with  a  known  voltage  applies  a  pressure  pulse  of 
known  amplitude  to  the  cavity  boundary,  propagating  compressional  waves  into  the 
surrounding  medium.  Lateral  heterogeneities  of  simple  geometries  (cylindrical  and  planar 
scattering  surfaces)  are  cast  into  the  surface  layer,  converting  incident  P-wave  energy  into 
shear  waves.  The  horizontally  polarized  shear  waves  (SH-waves)  trapped  in  the  surface 
layer  wave-guide  are  the  Love  Waves  we  will  measure  at  the  free-surface.  The  sensors  at 
the  surface  will  be  distributed  so  both  the  undisturbed  signal  and  the  signals  modified  by 
scattering  can  be  monitored  at  the  surface. 
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90  cm 


Figure  (a)  Model  with  planar  scattering 
surface. 


Figure  1 .  Schematic  of 
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Figure  (b)  Model  with  single  cylinder  or 
array  of  cylinders  as  scattering 
objects. 


-model  laboratory  experiment. 


PROGRESS 


An  analytic  solution  of  a  spherical  piezoelectric  crystal  v.cs  derived  to  predict 
expected  signal  levels  for  a  sensor  evaluation  experiment,  and  to  provide  a  boundary 
pressure  history  as  input  for  finite  element  code  calculation. 

SOLUTION  OF  A  SPHERICAL  PIEZOELECTRIC  SOURCE  IN  AN 
ELASTIC  MEDIUM 

The  source  is  a  spherical  shell  of  lead  zirconate  titanate  ceramic,  designated  PZT-4, 
of  outside  radius  a  =  0.6350  cm,  inside  radius  b  =  0.4826  cm,  and  thickness  h  =  0.1524 
cm.  The  electrical  polarization  is  radial  and  the  electrical  and  mechanical  properties  are 
spherically  symmetric  and  transversely  isotropic.  In  the  usual  spherical  coordinate  system 
(r,0,<}>),  the  properties  are  isotropic  with  respect  to  the  circumferential  coordinates  (0,<J>). 

The  source  is  cast  in  a  rock  simulant,  designated  Pourstone,  which  is  homogeneous 
and  isotropic.  Our  analysis  treats  the  source  subjected  to  a  voltage  pulse  and  an  interface 
reaction  from  the  sun-ounding  medium.  This  reaction  is  determined  by  an  analysis  of  the 
spherical  wave  propagation  in  the  surrounding  medium  caused  by  the  reaction  while 
ensuring  equal  radial  displacements  at  the  interface.  Approximations  are  introduced  to 
simplify  the  analysis  without  changing  the  basic  dynamic  response. 

PIEZOELECTRIC  SPHERICAL  SHELL 

The  standard  notations1  for  the  physical  quantities  involved  are 


Si 

strain 

Ti 

stress 

N/m2 

Em 

electric  field 

V/m  (N/C) 

SE 

elastic  compliance  coefficients 

m2/N 

CE 

elastic  stiffness  coefficients 

N/m2 

dmi 

piezoelectric  constants 

m/V  (C/N) 

Cmj 

piezoelectric  constants 

m/V  (C/N) 
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The  general  relationships  among  the  strain,  stress,  and  electric  field  are 


oi 


with 


Si  =  sfjTj  +  dniiEm 


(1) 


Tj  —  cjjSi  -  emjEm 
ij  =  1,2,.,., 6  m  =  1,2,3 


(2) 


Considerable  simplification  of  equations  (1)  and  (2)  results  when  we  make  use  of 
the  general  properties,  sy  =  Sjj  and  cy  =  Cjj,  transverse  isotropy  (rotational  symmetry  about 
the  X3  axis),  and  spherically  symmetric  electrical  excitation  with  Ej  =  E2  =  0.  Then, 


equations  (1)  and  (2)  become 

S2  *  Si  =  (sn  +  si2)Ti  +  S13T3  +  d3iE3  (3a) 

S3  =  2S13T1  +  S33T3  +  d33E3  (3b) 

T2  =  Ti  =  (cn  +ci2)Si  +C13S3  -  e3iE3  (4a) 

T3  =  2C13S1  +  C33S3  -  e33E3  (4b) 


The  local  (xj,  X2,  X3)  triad  is  the  spherical  coordinate  system  (8,0,r).  Thus,  stress 
components  Ti  and  T3  are  Tq  and  Tr,  strain  components  Si  and  S3  are  Sq  and  Sr,  and  the 
electric  field  component  E3  is  Ef. 

The  equation  of  motion  of  a  shell  element  is 


—  +  ^T3-T!)  =  p 


A 

at2 


(5) 


in  which  4(r»t)  is  the  outward  radial  displacement  and  p  is  the  density  of  the  PZT-4.  The 
strains  are 


Si=|  S3=^  (6) 

*  dr 

Substitution  in  (5)  of  the  stresses  from  (4)  followed  by  substitution  of  the  strains  from  (6) 
leads  to  the  governing  equation 
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(7) 


^2*  -\t  E 

C33  — —  +  2C33  7- - 2(Ci  1  +  C12  -  C13)  — 


dr 


dr 


3' l  dE  p 

P  — r  =  e33  —  +  2(e33  -  e3i)^ 
3t2  dr  r 


where  we  have  replaced  the  symbol  E3  by  E.  If  the  potential  at  the  inner  radius  r  =  b  is 
higher  than  the  potential  at  the  outer  radius  r  =  a  by  the  voltage  V,  the  field  is 


E(r,t)  = 

a  -  b 


m 

r2 


<K) 


Then  equation  (7)  becomes 


2 

2(cii +C12-C13)4-P“  =  -  e3iV(t)4 


dr 


dr 


dt 


a  -  b 


(9) 


At  the  free  inner  boundary,  r  =  b,  the  radial  stress  component  is  T3  =  0  so  by 
equations  (4b),  (6),  and  (8)  this  condition  is 


^  +  2-  =  522--ak-  •  ^ 

dr  r  C33  a  -  b  r2 


at  r  =  b 


(10) 


The  condition  at  the  outer  boundary  depends  on  the  problem  being  solved.  For  a  traction 
free  outer  boundary,  the  condition  is  T3  =  0  at  r  =  a,  that  is. 


t  2  j  -  e33  ab  .  V 
dr  r  C33a-b  r2 


at  r  =  a 


(ID 


To  obtain  an  upper  bound  on  the  interface  pressure  at  r  =  a,  we  have  Sj  =  S2  =  0  at  r  =  a, 
that  is, 

-  =  0  at  r  =  a 
r 

This  condition  corresponds  to  a  spherical  source  in  a  rigid  material.  Fc  our  case  of 
an  elastic  material,  the  condition  at  r  =  a  is  T3  =  -p(t),  the  interaction  pressure.  The  solution 
is  a  relationship  between  the  radial  displacement  at  r  =  a  and  the  interaction  pressure.  By 
analyzing  the  problem  of  a  pressure  p(t)  acting  in  a  spherical  cavity  of  radius  r  =  a  in  an 


5 


elastic  medium,  we  obtain  a  second  relationship  between  cavity  wall  radial  displacement 
and  the  ^rcscure  p(t).  Equating  the  displacements  gives  an  equation  for  the  required 
interface  pressure. 

The  governing  equation  (9)  and  the  various  boundary  conditions  can  be  solved 
explicitly  but  the  results  are  extremely  cumbersome.  Therefore,  for  aiding  the  interpretation 
of  experimental  results,  the  source  analysis  is  simplified  by  using  thin  shell  theory.  In  this 
theory,  the  radial  comp  nt  of  stress,  T3,  is  neglected. 

Let  the  mean  radius  of  the  shell  be  r  =  a  and  introduce  the  interface  pressure  p(t)  to 
be  determined.  The  shell  outward  radial  displacement  is  4(0  and  the  thickness  is  h.  The 
equation  of  motion  is 


P(t) 

h 


(13) 


and  if  in  equation  (3a)  we  set  T3  =  0  and  Si  =  4/a  to  provide  the  average  circumferential 
stress  Ti  in  equation  (13),  we  obtain 


4  +  co24  =  -  ^  +  C02ad3iE 


(14) 


where 

co2  =  2/pa2(sn  +  S12) 


(15) 


The  initial  conditions  of  interest  are 


4(0)  =  0  4(0)  =  0 

The  solution  of  equation  (14)  satisfying  the  initial  conditions  (16)  is 


co2ad3iE(T)  -  -±- 
ph 


sinco(t  -  x)dx 


If  the  driving  electric  field  is 


(16) 


(17) 


E(t)  =  E0(l  -e^“) 


(18) 
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the  shell  velocity,  according  to  the  displacement  (17),  is  given  by 


^  _  oc/co 


1  +  (a/co)2 


sincot  -  coscot  +  e  atj  - 


p(x)  cosco(t  -  t)dx 


(19) 


where 


v  =  coad3iE0 


(20) 


The  velocity,  v,  is  the  maximum  velocity  achieved  by  the  shell  in  a  vacuum  when  subjected 
to  a  step  voltage,  Eq.  Formula  (19)  gives  the  relationship  between  the  shell  radial  velocity 
and  the  interaction  pressure  from  the  surrounding  elastic  medium. 

SPHERICAL  CAVITY  IN  AN  ELASTIC  MEDIUM2 

When  the  stress-strain  relationships 


cr  =  (X  +  2ja)  —  +  2\  ~ 
dr  r 


oe  =  X  —  +  +  ^)  | 

dr 


are  substituted  in  the  equation  of  motion 


9ar  2  /  x  0  q 
—  +  f(ar-ae)  =  pm-f 

dr  dtz 


d\ 


we  obtain  the  displacement  equation 


2A=±.tA 

dr2  T  dr  T  c2  8t2 


(21) 


(22) 


(23) 


(24) 
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In  (21)  and  (22),  the  Or  and  Oe  are  the  radial  and  circumferential  stress  components,  and  X 
and  p.  are  the  Lame  constants  for  the  isotropic  elastic  medium.  In  (23),  pm  is  the  medium 
density.  In  (24),  c2  =  (X  +  2|i)/pm  where  c  is  the  elastic  wave  velocity. 

Introduction  of  the  displacer  ient  potential  defined  by 


reduces  (24)  to 


dV<P)  ..  i  d\rq) 

dr2  "c2  9t2 


The  solution  of  (26)  describing  outgoing  waves  is 

<p(r,0  =  }-  f(s)  s  =  t  -  (r  -  a)/c 


In  terms  of  the  function  f(s),  the  displacement,  velocity,  and  stresses  are 


£1  r 
at  'CT> 


in  which  f)  is  Poisson's  ratio. 

For  a  given  cavity  wall  velocity,  we  have,  from  (29),  the  equation 


f  (s) +  a  f  (s)  =  ca^'(s) 


(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 
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and  f(0)  =  f  (0)  =  0  at  the  wave  front  s  =  0.  The  solution  is 


f(s)  =  -  cae-05/®  |  e^a  ^(x)dx 


Jo 


with 


f(s)  =  - |  f(s)  -  ca£(s) 
f ’(s)  =  f(s)  -  ca^'(s)  +  c2£(s) 


(33) 


(34) 

(35) 


The  cavity  pressure  required  to  maintain  the  cavity  wall  velocity  is  found  by 
substituting  (33),  (34),  and  (35)  in  the  radial  stress  formula  (30)  and  setting  r  =  a  (s  =  t) 
and  Or(a,t)  =  -p(t).  These  steps  give 

p(t)  =  Pci<,)-^{f  +  5-M>5j 
with  f(t)  given  by  (33)  with  s  =  t. 

Formally,  the  velocity  of  the  interface  between  the  source  and  the  medium  is 
obtained  by  solving  the  equation  that  results  from  substituting  the  interface  pressure  (36)  in 
the  shell  velocity  equation  (19).  This  equation  is  inconvenient  so  we  introduce  an 
approximation  for  the  interface  pressure  (36).  By  integrating  f(t)  from  (33)  by  parts,  we 
obtain 


M  =  .  £  +  e(c/a)T 

a2 


I. 


e(c/a)T  £(t)dX 


(37) 


For  an  electric  field  excitation  (18)  on  the  sphere  in  a  vacuum,  the  second  term  on 
the  right  hand  side  of  (37)  is  much  smaller  than  the  first  term  if  a  «  0),  which  it  is  in  our 
case.  Assuming  a  similar  relative  magnitude  when  the  source  is  embedded  in  a  solid  allows 
us  to  approximate  f(t)  by 


a2 
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so  that  (36)  simplifies  to 


p(,)=pci«)+Hi«LJ2> 
a  1  --6 


4(0 


(38) 


Substitution  of  (38)  in  the  shell  equation  (14)  gives 

i  +  2yi  +  Cl\  =  cov(l  -  e-*) 

where 


(39) 


Y  =  Pn»c/2ph  (40) 

ft2  =  (02  +  (pmc2/pha)[2(l  -  2fi)/l  -  fi)]  (41) 

and  pm  is  the  density  of  the  medium.  The  wave  speed  in  the  medium  is  defined  by  c2  =  (k 
+  2pVpm.  In  terms  of  the  stiffness  coefficients,  the  angular  frequency  of  the  piezoelectric 
shell  in  a  vacuum  is  given  by 


to2  =  2(cu  +  C12  -  2c^3/C33)/pa2 

In  (39),  the  velocity,  v,  is  given  by  (20). 

The  solution  of  equation  (39)  satisfying  the  initial  conditions  (16)  is 


where 


°-S(0=  1  -  — 

tov  ’  T2 


Cl  +(y-a)2 


(Q2e-«'  -  ae-r 

(2y-a)cosQt  -  — — sinflt 

\ 

a 

(42) 


(43) 


n  -i1 


(44) 


For  excitations  with  rise  times  that  are  long  compared  to  the  natural  quarter  period, 
— 2 

we  have  a2  «  Q  and  the  radial  displacement  (43)  becomes 
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(45) 


•  ^(t)  =  1  -  c®1  -  ^  e-v  smQt 


The  shell  radial  velocity  obtained  by  differentiating  (45)  is. 


£(t)  =  e-01  -  e-if*  cosQt  -  =  sinQt 


Substitution  of  %(x)  from  (45)  into  (33)  and  performing  the  integration  gives 


r\2  . 

.  —  f(s)  =  1  -  —4 —  (£  e-“*  -  ae-<c/a)s) 

"2'"-  (c/a)  -  a  'a  ’ 


a(c/a) 


e-(c/a)s  .  e-YS  (cosQs  -  ^  sin£2s 


Formulas  (34)  and  (35)  determine  f  (s)  and  f'(s)  because  £(s),  £'(s),  and  f(s)  are  given  by 
(45),  (46),  and  (47).  Consequendy,  we  can  determine  £;,  £,  ar,  and  oe  by  (28)-(31). 


FLUID  MEDIUM 

If  we  set  ^  =  1/2  in  the  frequency  formula  (41), 

Q2  =  co2 

and,  by  (44) 


ft2  =  co2  -  Y2 


where  the  fluid  wave  velocity  is  c  =  (K/pm),  K  being  the  fluid  bulk  modulus.  In  our  case, 
Y  « (0  so  (45),  (46),  and  (47)  become 


^(s)  =  X-  ( 1  -  e-™  -  e'T*  sincos) 
to'  to  I 

4’(s)  =  |e*°“  -  e-75  (coscos  -  —  sincos JJ 
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-  f(s)  =  1  -  — } - (£e-«s  -  ae-W*)*} 

a2y  (c/a)  -  a  'a  ’ 

-  rcfofe)  |e-(c/i)s .  e-Y»  (coscos  -  —  sincosll 

0)2  '  '  (52) 

Stress  formulas  (30)  and  (31)  are  replaced  by  the  pressure  formula 

P  =  -  Pmf"(s)/r  (53) 

where  f’(s)  is  determined  by  (35),  (50),  (51),  and  (52). 

NUMERICAL  VALUES  (SOLID  MEDIUM) 

The  properties  we  require  of  the  PZT-4  ceramic  are 
Elastic  compliance  coefficients  (m2/N) 


cE 

all 

12.30  x  1012 

cE 

b12 

-  4.05  x  1012 

cE 

s13 

-  5.31  x  10-12 

oE 

o33 

15.50  x  10-12 

Elastic  stiffness  coefficients  (N/m2) 

rE 

cn 

13.90  x  1010 

C?2 

7.78  x  1010 

c?3 

7.43  x  10l° 

cf3 

1*1.50  x  1010 

Piezoelectric  constants  (m/V,  c/N) 

- 123  x  10-12 

d33 

289  x  10-12 

Piezoelectric  constants  (Nm/V,  c/m2) 

631 

-  5.20 

633 

15.10 

Density  p 

7.5  x  103  kg/m3 

The  properties  of  the  'Pourstone'  medium  are 

Young's  modulus  E  16.4  GPa 

Shear  modulus  H  6.3  GPa 
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Bulk  modulus 

K 

13.8  GPa 

Poisson's  ratio 

0 

0.3 

Lame  constant 

\ 

9.6  GPa 

Density 

Pm 

1.79  x  103  kg/m3 

P-wave  velocity 

°P 

3.52  mm/ps  (km/s) 

S-wave  velocity 

CS 

1.88  mm/(is  (km/s) 

nensions  of  the  spherical  source  are 

Outer  radius 

a 

6.350  mm 

Inner  radius 

b 

4.826  mm 

Thickness 

h 

1.524  mm 

The  driving  voltage,  V(t),  is  taken  in  the  form 

V(t)  =  V0(l-e-“t) 

applied  to  the  outside  of  the  spherical  shell.  If  the  value  at  the  midradius,  (a  +  b)/2,  is 
chosen  to  represent  the  field  strength,  then  according  to  (8)  and  (54), 

E(t)  =  -  644  Vo(l  -  e-01)  volts/m 

By  letting  the  voltage  reach  90%  of  V0  in  20  (is,  the  value  of  a  is  determined  as  a  =  0.1 15 
(is’1.  Also,  if  V0  =  300  volts,  we  have  E0  =  -0. 1932  x  106  volts/m.  Hence, 

E(t)  =  -  0.1932  x  106(1  -e-°-115t) 

where  t  has  |is  units. 

The  natural  angular  frequency  of  the  spherical  source  (in  a  vacuum),  according  to 
(15),  is  0)  =  0.895  x  106  rad/s  and  the  natural  frequency  is  f  =  co/2jt  =  142  kHz. 

If  the  voltage  is  applied  slowly  to  the  free  shell,  the  maximum  radial  displacement  is 
=  ad3iE0  =  0.15  x  10'3  mm  =  0.15  (im  . 

If  the  same  voltage  (300  volts)  is  applied  instantaneous,  the  maximum  displacement  is  2£s, 
and  the  maximum  velocity  is 

max  £  =  v  =  coad3iE0  =  (o£s  =  13.4  cm/s 


(54) 
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If  the  voltage  is  applied  slowly  to  the  shell  confined  by  a  rigid  medium,  the  interface 
pressure  is  p  =  -  (2h/a)Ti  where  the  circumferential  stress  is  Tj  =  -  d3iEo/(sn  +  S12X  The 
magnitude  of  this  pressure  is  p  =  14.05  bars  (206  psi). 

In  Equation  (39)  governing  theshell  motion  in  an  elastic  medium,  the  numerical 
values  of  y  and  £2,  and  consequently  &  are 

y  =  0.276  x  106  rad/s  £2  =  1.073  x  106  rad/s  £2  =  1.037  x  106  rad/s 

— 2 

For  comparison,  we  note  that  a  =  0.1 15  x  106  s_1.  Because  a2  «  £2  formulas  (45), 

(46),  and  (47)  for  the  shell  radial  displacement  and  velocity  and  the  potential  function  f(s) 
are  applicable.  The  value  of  c/a  occurring  in  f(s)  is  c/a  =  0.554  x  106  s*1. 

The  oscillatory  part  of  the  solution  has  a  frequency  of  f  =  £2/2 71  =  165  kHz  with  a 
period  of  T  =  1/f  =  6  ps. 


NUMERICAL  VALUES  (WATER) 

Using  a  density  of  pm  =  1  gram/cm3 
kbar),  we  obtain 

Wave  velocity 
Damping  parameter 
Spring  parameter 
Angular  frequency 
Velocity  parameter 
Excitation  parameter 


and  a  bulk  modulus  of  K  =  2.245  GPa  (22.45 

c  =  (K/pm)1/2  =  1.50  mm/ps 
y  =  0.066  x  1 06  rad/s 
£2  =  0)  =  0.895  x  106  rad/s 
£2  =  0.893  x  106  rad/s®  £2 
c/a  =  0.236  x  106  s'1 
a  =  0.115  x  106  s1 


APPLICATION  OF  SOLUTION  IN  A  WATER  MEDIUM 

We  applied  the  solution  just  derived  for  the  case  of  a  spherical  piezoelectric  source 
in  a  water  medium,  and  compared  the  calculated  and  measured  pressure  histories  at  three 
radii  from  the  source.  The  experimental  configuration  is  shown  in  Figure  2.  In  the 
experiments,  the  crystal  was  excited  with  a  known  voltage  history,  and  free-field  water 
pressure  histories  were  measured  at  radii  of  0.91, 1.51,  and  3.38-cm  measured  from  the 
center  of  the  source.  The  input  voltage  history  to  the  crystal  is  shown  in  Figure  3.  The 
voltage  reaches  90  %  of  the  peak  value  of  316  volts  at  10  ps,  so  in  the  formulation 
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RA  M  2336-2A 


Figure  2.  Configuration  for  measuring  stress  pulse  amplitudes  in  water  at  different 
radii  from  the  source  and  sphericity  of  piezoelectric  source. 


5 


Figure  3.  Input  voltage  history  to  the  piezoelectric  source  in  the  water 
pressure  experiment. 
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described  previously,  the  value  of  a  is  0.230  x  106  (is*1,  and  the  field  strength  as  a 
function  of  time  is  described  by: 

E(t)  =  -0.2075  x  106(1  -e-°-23  *) 

where  the  field  strength,  E(t),  has  units  Volts/m  and  t  has  units  |is. 

The  calculated  and  measured  pressure  histories  are  shown  superposed  for  the  three 
gage  locations  in  Figure  4,  5,  and  6.  At  the  first  gage  location,  the  measured  and  calculated 
histories  show  very  good  agreement  in  oscillation  frequency  and  peak  amplitude.  At  the 
further  out  locations,  the  measurements  show  larger  damping  than  predicted  by  the 
calculation,  but  satisfactory  agreement  is  observed  in  peak  pressure  and  the  oscillation 
frequency. 

APPLICATION  OF  SOLUTION  TO  AN  ELASTIC  MEDIUM 

The  agreement  between  calculation  and  experiment  in  a  water  medium  is  sufficient 
to  extend  the  solution  to  an  elastic  medium  and  estimate  the  expected  amplitudes  of  velocity 
and  displacement  at  different  ranges  in  the  medium.  In  addition,  the  calculated  pressure 
history  at  the  source/medium  interface  can  be  used  as  an  input  boundary  condition  for  finite 
element  calculations  of  the  experiment  shown  in  Figure  1.  The  results  of  the  finite  element 
calculations  will  be  useful  for  estimating  expected  signal  levels  on  the  surface  and  assist  is 
instrumentation  selection  for  that  experiment.  The  elastic  medium  is  a  rock  simulant  called 
pourstone,  with  the  elastic  properties  listed  on  page  (16).  Using  the  same  voltage  history 
shown  in  Figure  3  as  the  input,  the  calculated  velocity  histories  at  radii  of  0.635  cm 
(ujurce/medium  interface),  1.5-  and  2.5-cm  are  shown  superposed  in  Figure  7.  The 
corresponding  displacements  for  the  three  locations  are  shown  in  Figure  8,  and  the  radial 
and  circumferential  stress  histories  are  shown  in  Figures  9  and  10,  respectively. 

EVALUATION  EXPERIMENT  IN  POURSTONE 

We  performed  experiments  to  compare  the  calculated  and  measured  velocity  and 
displacement  histories  in  a  sphere  of  pourstone.  The  experimental  configuration  is  shown 
in  Figure  1 1.  In  this  experiment,  copper  loops  were  cast  in  the  midplane  of  a  12-cm 
diameter  sphere  of  pourstone  to  measure  radial  particle  velocity  histories  at  different  radii 
from  the  source.  The  specimen  is  placed  in  an  external  magnetic  field,  and  we  measure  the 
induced  voltage  as  the  conductor  cuts  flux  lines  during  passage  of  the  stress  wave.  The 
particle  velocity  is  proportional  to  the  induced  voltage,  the  magnetic  field  strength,  and  the 
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Figure  4.  Comparison  of  measured  and  calculated  pressure  histories 
in  water  at  a  range  of  0.91  -cm  from  the  center  of  the  source. 
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RA  7206  5 


Figure  5  Comparison  of  measured  and  calculated  pressure  histories 
in  water  at  a  range  of  1 .51 -cm  from  the  center  of  the  source. 
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Figure  6.  Comparison  of  measured  and  calculated  pressure  histories 
in  water  at  a  range  of  3.38-cm  from  the  center  of  the  source 
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Figure  7.  Calculated  velocity  histories  at  3  ranges  in  pourstone  from  a 
spherical  piezoelectric  crystal  excited  by  316  volts. 
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Figure  8.  Calculated  displacement  histories  at  three  ranges  in  pourstone  from 
a  spherical  piezoelectric  crystal  excited  by  316  volts. 
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Figure  9.  Calculated  radial  stress  histories  at  three  ranges  in  pourstone  from  a 
spherical  piezoelectric  crystal  excited  by  316  volts. 
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RA- 7206- 10 

Figure  1 0.  Calculated  circumferential  stress  histories  at  three  ranges  in 

pourstone  from  a  spherical  piezoelectric  crystal  excited  by  316  volts. 
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*  Particle  Velocity  Measurements  Only 


Figure  1 1 .  Configuration  for  source/sensor  evaluation  experiment. 


length  of  the  conductor.  To  increase  the  expected  signal,  we  increased  the  gage  length  by 
constructing  gages  consisting  of  10  windings  of  copper  wire.  Accelerometers  and  PVDF 
foils  were  mounted  on  the  surface  of  the  sphere  to  measure  free-surface  acceleration  and 
strain,  respectively. 

We  are  currently  analyzing  the  data  obtained  from  the  particle  velocity  gages.  The 
accelerometers  and  PVDF  foils  did  not  produce  satisfactory  signals  because  the  range  of  the 
measurement  was  too  large.  Consequently,  we  plan  to  fabricate  a  smaller  specimen  to 
evaluate  these  gages  at  locations  closer  to  the  source. 
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